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Abstract We study the cooling of hybrid stars coupling with spin-down. Due to the spin- 
down of hybrid stars, the interior density continuously increases, different neutrino reac- 
tions may be triggered(from the modified Urea process to the quark and nucleon direct 
Urea process) at different stages of evolution. We calculate the rate of neutrino emissivity 
of different reactions and simulate the cooling curves of the rotational hybrid stars. The 
results show the cooling curves of hybrid stars clearly depend on magnetic field if the 
direct urea reactions occur during the spin-down. Comparing the results of the rotational 
star model with the transitional static model, we find the cooling behavior of rotational 
model is more complicated, the temperature of star is higher, especially when direct urea 
reactions appear in process of rotation. And then we find that the predicted temperatures 
of some rotating hybrid stars are compatible with the pulsar's data which are contradiction 
with the results of transitional method. 
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1 INTRODUCTION 

The interior of neutron stars contain matter beyond the nuclear saturation density which we have not 
seen yet. The cooling of neutron stars give us an important tool to study the properties of such dense 
matter. The traditional investigations of cooling often adopt the static star model which is not connect 
with spin-down. It is well known that neutron stars would spin-down due to magnetic dipole radiation. 
The spin-down compression of stars may lead to the changes in chemical composition(from nucleon 
matter to deconfined quark matter) and structure(mixed phase and quark phase appearing). The coupling 
of cooling and spin-down correlates stellar surface temperatures with rotational state as well as time. 
The interlinked processes of spin-down and cooling present intriguing prospects of gain insight in the 
fundamental properties of dens e mat ter in neutron stars by confrontation with thermal emission data 
from observations(Stejner et al ( l2008b ). (Page et al. (l2006l) .Yakovlev & Pethick (|2004 ). 

Neutron stars are born with temperatures above The dominant cooling mechanism of stars 

is the neutrino emission from the interior for the first several thousand years after birth, which can be 
generated via numerous reactions(Yakovlev & Pethick (|2004|) ). For the nucleon direct urca(NDU) re- 
action, the most efficient one, is only possible if the fraction of proton exceeds a certain threshold. It 
is impossible to satisfy conservation of momentum unless the proton fraction exce ed the value where 
both charge neutrality and the triangle inequality can be observed (Lattimer et al. (ll99lb ). Hence the 
traditional investigations of neutron stars often divide the cooling process into two regimes that are slow 
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and fast cooling due to slow and fast neutrino emission respectively. Slow cooling occurs in the low 
mass stars via neutrino emission produced mainly by the nucleon modified Urca(NMU) process. The 
fast cooling occurs in stars with mass critical one Moiii is of course model dependent) via the NDU 
process. Comparing with the observed data, we can see that the fast cooling process would result in con- 
tradiction between the predicted te mpera tures of stars and the observations, which gives the challenge 
to the traditional model. (Page et al. 

The neutron stars containing quark matter are called hybrid stars. Hybrid stars have more compli- 
cated interior structure and matter composition than the purely neutron stars. Due to spin-down com- 
pression, the interior density gradually increases, some the thermodynamic quantities such as neutrino 
emission luminosity and total heat capacity of stars continuously change with rotational frequency. 
Especially, the rate of neutrino emission would have an abrupt rise because of direct urea reactions oc- 
curring which would induce the rapid fall of temperatures of hybrid stars. Appearing of different direct 
Urea processes can result in different cooling behavior in the stages of evolution of the stars. The main 
difference between our model and traditional model is that our model combine the equation of thermal 
balance with the rotational stars structure and take magnetic dipole radiation model to investigate the 
changes of thermodynamic quantities and the cooling behavior of hybrid stars with spin frequency as 
well as time. The simulation of cooling curves of rotational stars are more complex than the traditional 
cases because of the magnetic field dependence and changes of rotational state. Surface temperatures of 
some stars including fast cooling processes are compatible with the pulsar's data due to spin-down. 

We take Glendenning's hybrid stars model (Glendenning (1 1997b ) based on the perturbation 
theory(Hartle(jT967),Chubarian et al (2000)) to study the rotational structure of stars. In our calcula- 
tion, we only choose the simplest nucleon matter composition, namely neutrons, protons, electrons, and 
muons, and ignore superfluidity and superconductivity. 



2 HYBRID STARS 



As the stars spin-down, the nuclear matter are continuously converted into quark matter by the exother- 
mic reactions, i.e. n ^ u + 2d, p 2u + d, s quarks immediately appear after weak decay. The decon- 
finement phase transition brought forward by Glendenning ( (J 997) . 11992)) who had realized firstly the 
possibility of the occurrence of a mixed phase(MP) of hadron matter and quark matter in a finite den- 
sity range inside neutron stars. Global charge neutrality of MP can be achieved by a positively charged 
amount of hadronic matter and a negatively charged amount of quark matter. The Gibbs condition for 
mechanical and chemical equilibrium at zero temperature between the two phases reads 

PHpif^n, IJ^e) = PQpif^n, 1-l.e) ^PMP- (1) 

here /i„ is the chemical potential of neutron and /ig is the chemical potential of the electron. The condi- 
tion of global charge neutrality in the MP is 

Q 

= — = XIQP + (1 - XJQHP- (2) 

Here x = ^q/^ is the quark fraction in the MP, thus the energy density emp of the MP follows as 

E 

ej\/p = y = X^QP + (1 - X)^HP- (3) 

We can obtain the equation of state (EOS) of MP using above equations. For the hadron part of star, 
we adopt the Argonne Vli + 5v + UIX* model (Akmal et al. ( Il998h ) which is based on the models 
for the nucleon interaction with the inclusion of a parameterized three-body force and relativistic boost 
correc tions. For the quark matter, we use the EOS of a effective mass bag-model EOS(Schertle et al. 
(Il997h ). In Figure 1, we show the EOS for hybrid stars with deconfinement transition as described 
above. We choose the parameters for quark matter EOS with s quark mass = IhOMeV, coupling 
constant g = 'i and different bag constants B = 108Af eVfm~^ and B = 136MeVfm~^. 
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Fig. 1 Model EOS for the pressure of hybrid star matter as a function of the energy density. 
The hadronic phase EOS is Argonne 1^18 + Sv + U IX* model. The quark phase is effective 
mass MIT bag model with bag constant B=136 and 108 MeV fm^^ respectively. 

The rotati on of stars lead to the change in structure. In present work, we apply Hartle's approach 
(Hartle (Il967h ') to investigate the rotational structure of the stars. It is based on the treatment of a ro- 
tating star as a perturbation of non-rotating star, which can be obtained by expanding the metric of an 
axially symmetric rotating star in even powers of the angular velocity Vl. We assu me th e frequency of 
stars a t birth close to the Keple r(mass-shedding limit) frequency(Hartle & thorne ( ll968h .Lattimer et al 
( I1990I) . De Araujo et al d 1995b . We have shown the central density of rotating stars of different grav- 
itational mass, as a function of its rotational frequency in Figure 2. These sequences of star models 
have same total baryon number but different central density and angular velocity. It is well known the 
perturbative approach fails when the angular velocity approaches the mass-shedding limit. However, the 
rotational frequencies of all known pulsars are much lower than the mass-shedding limit. So we can use 
the perturbative theory to investigate the structure of rotating stars unless we are specifically interested 
in the mass-shedding problem(Benhar et al (l2005i) ). In Figure 2, dotted horizontal lines indicate decon- 
fined quark matter produced and dashed horizontal lines indicate NDU nucleon direct Urea processes 
triggered. The appearance of quark matter will result in the occurrence the quark direct Urca(QDU) 
processes. Such processes are also efficient, but somewhat weaker than the NDU ones. We observe 
spin-down of stars lead to the changes in chemical composition and structure in Figure 2. In the case 
of bag constant B = 108A/ey/TO^'^(left penal), the quark matter is produced in the interior for the 
stars mass M > 1.4il/0(here denote the static mass) in process of rotation. For M = I.6M0 star, for 
example, a purely neutron star can be transformed into a hybrid star and NDU reactions are triggered 
during spin-down. It is known that the occurrence of QDU and NDU processes lead to the rapid fall 
of temperatures of stars. In the section four, we discuss in detail spin-down leads to changes of the 
temperatures and cooling curves. 

3 NEUTRINO EMISSIVITIES 

The cooling of stellar could be carried out via two channels - neutrinos emission from the entire star 
and thermal emission of photons through transport of heat from the internal layers to the surface. The 
emission of neutrinos could carry away energy which would provide efficient cooling for warm neutron 
stars. 
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Fig. 2 Central density as a function of rotation frequency for rotating hybrid stars with dif- 
ferent gravitational mass at zero frequency. All sequences have constant total baryon number. 
Dotted horizontal lines indicate deconfined quark matter is produced and dashed horizontal 
lines indicate nucleon direct Urea process is triggered. 

For the hadron matter, we mainly consider three kinds of processes which are NDU, NMU as well 
as nucleon bremsstrahlung(NB). For the quark matter, we take QDU processes of unpaired quarks, the 
quark modified Urea (QMU) processes and the quark bremsstrahlung(QB) are considered. The most 
powerful neutrino emission is provided by direct Urea process. The rate of neutrino emissivity of NDU 
processes n pev, pe nv is given by 

ENDU - 4.0 X lO'^Y.^f/^TlQt ergcm-h-^ (4) 

Ps 

where Tg is the temperature in units of 10^ K, ps — 0.16/m~^ is the nuclear saturation density, Qt — 
Oj PFe. + PFr, — PFn) IS the threshold factor, with 9{x) being 1 for .t > and zero otherwise(Lattimer et 
al. (ll99ll) ). NDU processes can occur when the fraction of proton exceeds 11%. The presence of muons 
would raise it to about 15%. 

The neutrino emissivity of quark matter has been given firsdy by IwamotodHH). The QDU pro- 
cesses have been estimated as 

EQDU - 8.8 X lO'^ac (^y^ Y}/^T^ erg cm'^ s"! (5) 

with the standard value of the QCD coupling constant ac — 0.1. 

The emissivity of NMU and NB processes in the non-superfluid npe matter are usually taken from 
Friman& Maxwell d 19791) . in which the one-pion-exchange Born approximation with phenomenological 
corrections was us ed for consideration. For the emissivity of QMU and QB processes, we take the results 
oflwamotodHH. 

4 COOLING CURVES WITH SPIN DOWN 

The traditional standard cooling model which often based on the Tolman- Oppen heimer-Volkoff(TOV) 
equation of hydrostatic equilibrium(Page et al. (l2006l) .Yakovlev & Pethick (|2004|) ). All thermodynamic 
quantities (neutrino emission luminosity and total heat capacity etc) are calculated when rotational fre- 
quency of stars is zero. 
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In present work, we investigate the cooling of hybrid stars with spin-down. The equation of ther- 
mal balance has been assumed spherical symmetry although it may be broken in a rotating star It is 
reasonable for slowly rotating stars wh ich co uld be treated as a perturbation to change the structure 
and chemical composition(Stejner et al ( |2008|) ). We combine the equati on of thermal balance with the 
rotating structure equations of the stars(kang & Zheng (l2007b . Hartle d 19671) ) and rewrite the energy 
equation in the approximation of isothermal interior(Glen & Sutherland (Il980h ) 



Cv{T.,v)^ = -L^{T,,v)~L?;{Ts,v) (6) 



L^iT,, v) = / e{r, r)(l - ±^)~^/^e^'^A7:r^dr (8) 







Where Ts is the effective surface temperature, Ti{t) ~ T{r,t) is the red-shifted internal tempera- 
ture; T{r,t) is the local internal temperat ure of matter, and $(r) is the metric function(describing 
gravitational red-shift)(Yakovlev & Haesel ( l2003b ). Furthermore L'^{Ti, v)and Cv{Ti,v) are the to- 
tal red-shifted neutrino luminosity and the total stellar heat capacity respectively which are functions 
of rotation frequency v and temperature T; c(r, T) is the heat capacity per unit volume. = 
ATiR^{v)(TTg{l — Rg/R) is the surface photon luminosity as detected by a distant observer(i?g is 
the gravitational radius of stellar ). The effective surface temperature detected by a distant observer 
is = T,^\/\ ~ Rn / R. Ts is obtained from the i nternal temperature by assuming an envelope 
model(Gudmundsson et al. (ll983h .Potekhin et al (ll997h ').The spin-down of stars is due to the magnetic 
dipole radiation. The evolution of rotation frequency is given by 

where / is the stellar moment of inertia, fi — ^B„iR^ is the magnetic dipole moment, and 6 is the 
inclination angle between magnetic and rotational axes. According to the Eq.(6)-(9), we can simulate 
the cooling of hybrid stars during spin-down. 

In Figure. 3 we present the cooling curves of a \.61\Iq rotational hybrid star for different magnetic 
field strengths (10^ — 10^'^G) with bag constant B = \Q^MeV fm~^ . We find the cooling curves of 
the star are different in different magnetic fields, especially for the strong magnetic fields cases. As 
we know, the stronger magnetic field is, the faster rotational frequency slows down due to the magnetic 
dipole radiation, the direct Urea processes are triggered at the earlier time, which result in rapidly cooling 
in a shorter time. In the cases of weaker fields(i? < 10^°G), direct urea processes appear about at the 
stage of photons cooling, the effect of magnetic fields is not important at the era. Using the same model 
in Figure 3, we show the neutrino emissivity of different reactions as well as the photon luminosity as 
functions of time and rotational frequency with magnetic field i?,„ = lO^^G in Figure 4. As spin-down 
of the star, deconfined quark matter appear in the core of star and then QDU reactions are triggered at 
spin frequency v = 1 123Hz. From then on, the neutrino emissivity of direct urea processes dominate the 
cooling of star. For star in age 10^'^ < t < lO^^yr, the QDU reactions provide most efficient neutrino 
emissivity. We find NDU processes are triggered at spin frequency v ~ 492Hz which lead to the rapid 
increase of neutrino emissivity. The photon emission control the cooling curves from about t = lO^yr 
after birth. 

In Figure 5, we plot the cooling curves of hybrid stars for different mass. The results of the ro- 
tational model with magnetic field B,n = 10^ and the traditional model are both presented in the 
picture(The per star model include two cooling curves.). The results show the temperatures of the ro- 
tating stars model are higher than the transitional static cases, especially for the star including NDU 
reactions(1.6 and 1.7 Mq stars for bag constant B = 108MeVfm~^, 1.8 and 1.9 Mq stars for bag 
constant B = 136M eVfm~^). As we know that the traditional model use the equation of hydrostatic 
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Fig. 3 Cooling curves of 1.6 Mq rotational hybrid star for various magnetic fields (curve a: 
IQi^G, b: IQi^G, c: lO^^G, d: lO^^G, e: lO^G). 
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Fig. 4 1.6 Mq hybrid star with magnetic field B„i = 10^^ for bag constant equal 108 
MeVfm^^, the neutrino emissivity of different reactions as well as the photon luminosity as 
functions of time(left panel) and rotation frequency(right panel). 

equilibrium to study the configuration of stars. The central density of stars and the thermodynamic 
quantities of neutrino emission luminosity and total heat capacity etc are calculated with rotational fre- 
quency of stars being zero. For static model, these hybrid stars are born when direct reactions occur. 
For rotational model, direct urea processes appear with the central density gradually increases during 
spin-down of stars, which induce the temperature of rotating model is higher than the static model at 
the stage of neutrino cooling. From Figure 5, we can observe the cooling curves of some rotating hy- 
brid stars(1.6MQand 1.7Mq in left panel, 1.8MQandl.9MQ in right panel) are compatible with the 
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Fig. 5 Cooling curves of hybrid stars for different mass with magnetic field = lO^^G and 
the cases in absence of rotation. The curves of the higher temperature correspond to our model 
for different mass stars. The observational data 1 to 11 are taken from tables 1 and 2 of Page 
et al.t2004). These stars are; 1. RX J0822-4247, 2. IE 1207.4-5209, 3. PSR 0538+2817, 4. 
RX J0002+6246, 5. PSR 1706-44, 6. PSR 0833-45, 7. PSR 1055-52, 8. PSR 0656+14, 9. PSR 
0633+1748, 10. RX J1856.5-3754, and 11. RX J0720.4-3125.. The next two stars, labeled as 
a and b: a. PSR B0531+21(Weisskopf et al.([2004))), b. PSR J0205+6449(Slane et al.(|2003)) 



observed data(Page et al. (l2004 . Weisskopf et al. (l2004 ,Slane et al. which are contradiction with 

the cooling curves of the transitional model. We find hybrid stars have complicated cooling behavior 
during spin-down. For M = 1.5A/0 hybrid star(left panel), the effect of spin-down is weaker than the 
cases of M = 1.6Mq and I.TMq because the interior of M = 1.5Mq star only appear QDU reactions, 
but M = 1.6Mq and 1.7Mq stars include QDU and NDU reactions in process of rotation. Comparing 
M = I.6M0 with I.TMq stars(left panel), we find earlier appearing of NDU processes lead to more 
rapid cooling for the rotating hybrid stars. 

We take the value of bag constant of relatively large(i3 = 108, 136MeVfm~^) in the paper. 
Comparing the results for two parameters, we find that critical mass of hybrid stars is smaller, the effect 
of spin-down is more important for the smaller mass hybrid stars with the decrease of bag constant. 
In the case of the smaller bag constant, spin-down may lead to the larger changes for temperature of 
the low mass hybrid stars. However, cooling behavior of hybrid stars for the different bag constant are 
similar if the changes of the rotational state are similar during spin-down process. 

5 CONCLUSIONS AND DISCUSSIONS 

The cooling of hybrid stars with spin-down have been studied in the paper. We combine the cooling 
equation with the rotational structure equation of stars and simulate the cooling curves of the hybrid 
stars. The results show the cooling curves have a clear magnetic field dependence if QDU or/and NDU 
reactions are triggered during spin-down. The time of direct Urea reactions triggered controls the occur- 
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ring of fast cooling. Comparing the cooling curves of the rotational models with the static models, we 
find the cooling behavior of rotational models are more complicated, the temperature of stars are higher, 
especially when direct urea reactions appear in process of rotation. We also find that the cooling curves 
of some rotational hybrid stars are consistent with the observational data in case of direct Urea neutrino 
emission. Through considering the inclusions of superfluidity and superconductivity, we expect cooling 
curves of these rotational hybrid stars can match the pulsar data well in future study. 

It has been studied by many investigators that different EOS of hadron phase and model parameters 
of quark phase (bag constant B, coupling constant g) would influence the phase transitio n den sities, 
rotational structure of hybrid stars and corresponding internal structure etc(Schertler et al. (l2000h : Pan 
et al.(2006). We will investigate the effect of these parameters on the cooling of rotating hybrid stars in 
detail. 

This work is supported by NFSC under Grant Nos. 10747126. 
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